New Approach for Temporal Stability Evaluation of Pseudo-Invariant Calibration Sites (PICS)

Pseudo-Invariant Calibration Sites (PICS) are one of the most popular methods for in-flight vicarious radiometric calibration of Earth remote sensing satellites. The fundamental question of PICS temporal stability has not been adequately addressed. However, the main purpose of this work is to evaluate the temporal stability of a few PICS using a new approach. The analysis was performed over six PICS (Libya 1, Libya 4, Niger 1, Niger 2, Egypt 1 and Sudan 1). The concept of a Virtual Constellation was developed to provide greater temporal coverage and also to overcome the dependence limitation of any specific characteristic derived from one particular sensor. TOA reflectance data from four sensors consistently demonstrating stable calibration to within 5%the Landsat 7 ETM+ (Enhanced Thematic Mapper Plus), Landsat 8 OLI (Operational Land Imager), Terra MODIS (Moderate Resolution Imaging Spectroradiometer) and Sentinel-2A MSI (Multispectral Instrument)were merged into a seamless dataset. Instead of using the traditional method of trend analysis (Students T test), a nonparametric Seasonal Mann-Kendall test was used for determining the PICS stability. The analysis results indicate that Libya 4 and Egypt 1 do not exhibit any monotonic trend in six reflective solar bands common to all of the studied sensors, indicating temporal stability. A decreasing monotonic trend was statistically detected in all bands, except SWIR 2, for Sudan 1 and the Green and Red bands for Niger 1. An increasing trend was detected in the Blue band for Niger 2 and the NIR band for Libya 1. These results do not suggest abandoning PICS as a viable calibration source. Rather, they indicate that PICS temporal stability cannot be assumed and should be regularly monitored as part of the sensor calibration process

Evaluation of Low-Cost Radiometer for Surface Reflectance Retrieval and Orbital Sensor’s Validation

This paper evaluates the Arable Mark 2 sensor, an automated and low-cost radiometer, for its potential to retrieve surface reflectance data and validate orbital sensors such as the Landsat-8 (L8) Operational Land Imager (OLI) Level 2 product. While orbital sensors are widely used for monitoring solar radiation changes, managing natural resources, and understanding climatic trends, atmospheric effects can make it challenging to obtain accurate measurements. Equipped with multiple sensors, including long-wave and short-wave radiometers, the Arable Mark 2 sensor can measure upwelling and downwelling irradiance to calculate surface reflectance. To assess the accuracy and consistency of the Arable Mark 2 sensor, the study performed a cross-calibration using a ground truth measurement collected with the Analytical Spectral Device (ASD) as the reference point. Additionally, a spectral band adjustment factor (SBAF) was applied across the calibrated Arable surface reflectance to compare it against the orbital sensor. An automated library aided in calculating SBAF for the days with unavailable hyperspectral data. The study found that the Arable Mark 2 sensor can provide accurate surface reflectance data that can be used for orbital sensor validation. The Arable sensor was successfully calibrated against the ASD FieldSpec with an average difference of less than 1/10 reflectance unit (reflectance unit = 0.01) for the blue, green, yellow, and red bands. The red-edge and NIR-1 bands showed an average difference of less than 1/2 reflectance units, while the NIR-2 band had an average difference of less than 1/10 reflectance unit of calibration accuracy. The calibrated Arable surface reflectance data was then compared against orbital sensor surface reflectance data, and the results showed good agreement between the two datasets. The study concludes that the low-cost and automated nature of the Arable Mark 2 sensor makes it a promising tool for surface reflectance retrieval and orbital sensor validation

Extended Pseudo Invariant Calibration Site-Based Trend-to-Trend Cross-Calibration of Optical Satellite Sensors

Satellite sensors have been extremely useful and are in massive demand in the understanding of the Earth’s surface and monitoring of changes. For quantitative analysis and acquiring consistent measurements, absolute radiometric calibration is necessary. The most common vicarious approach of radiometric calibration is cross-calibration, which helps to tie all the sensors to a common radiometric scale for consistent measurement. One of the traditional methods of cross-calibration is performed using temporally and spectrally stable pseudo-invariant calibration sites (PICS). This technique is limited by adequate cloud-free acquisitions for cross-calibration which would require a longer time to study the differences in sensor measurements. To address the limitation of traditional PICS-based approaches and to increase the cross-calibration opportunity for quickly achieving high-quality results, the approach presented here is based on using extended pseudo invariant calibration sites (EPICS) over North Africa. With the EPICS-based approach, the area of extent of the cross-calibration site covers a large portion of the North African continent. With targets this large, many sensors should image some portion of EPICS nearlydaily, allowing for evaluation of performance with much greater frequency. By using these near-daily measurements, trends of the sensor’s performance are then used to evaluate sensor-to-sensor daily cross-calibration. With the use of the proposed methodology, the dataset for cross-calibration is increased by an order of magnitude compared to traditional approaches, resulting in the differences between any two sensors being detected within a much shorter time. Using this new trend in trend cross-calibration approaches, gains were evaluated for Landsat 7/8 and Sentinel 2A/B, with the results showing that the sensors are calibrated within 2.5% (within less than 8% uncertainty) or better for all sensor pairs, which is consistent with what the traditional PICS-based approach detects. The proposed cross-calibration technique is useful to cross-calibrate any two sensors without the requirement of any coincident or near-coincident scene pairs, while still achieving results similar to traditional approaches in a short time

Evaluation of RadCalNet Output Data Using Landsat 7, Landsat 8, Sentinel 2A, and Sentinel 2B Sensors

In 2013, the Committee on Earth Observation Satellites (CEOS) Working Group on Calibration and Validation (WGCV) Infrared and Visible Optical Sensors Subgroup (IVOS) established the Radiometric Calibration Network (RadCalNet), consisting of four international test sites providing automated in situ measurements and estimates of propagated top-of-atmosphere (TOA) reflectance. This work evaluates the ‘reliability’ of RadCalNet TOA reflectance data at three of these sites—RVUS, LCFR, and GONA—using Landsat 7 ETM+, Landsat 8 operational land imager (OLI), and Sentinel 2A/2B (S2A/S2B) MSI TOA reflectance data. This work identified a viewing angle effect in the MSI data at the RVUS and LCFR sites; when corrected, the overall standard deviation in relative reflectance differences decreased by approximately 2% and 0.5% at the RVUS and LCFR sites, respectively. Overall, the relative mean differences between the RadCalNet surface data and sensor data for the RVUS and GONA sites are within 5% for ETM+, OLI, and S2A MSI, with an approximately 2% higher difference in the S2B MSI data at the RVUS site. The LCFR site is different from the other two sites, with relative mean differences ranging from approximately -10% to 1%, even after performing the viewing angle effect correction on the MSI data. The data from RadCalNet are easy to acquire and use. More effort is needed to better understand the behavior at LCFR. One significant improvement on the accuracy of the RadCalNet data might be the development of a site-specific BRDF characterization and correction

Validation of Expanded Trend-to-Trend Cross-Calibration Technique and Its Application to Global Scale

The expanded Trend-to-Trend (T2T) cross-calibration technique has the potential to calibrate two sensors in much less time and provides trends on a daily assessment basis. The trend obtained from the expanded technique aids in evaluating the differences between satellite sensors. Therefore, this technique was validated with several trusted cross-calibration techniques to evaluate its accuracy. Initially, the expanded T2T technique was validated with three independent RadCaTS RRV, DIMITRI-PICS, and APICS models, and results show a 1% average difference with other models over all bands. Further, this technique was validated with other SDSU techniques to calibrate the newly launched satellite Landsat 9 with 8, demonstrating good agreement in all bands within 0.5%. This technique was also validated for Terra MODIS and ETM+, showing consistency within 1% for all bands compared to four PICS sites. Additionally, the T2T technique was applied to a global scale using EPICS Global sites. The expanded T2T cross-calibration gain result obtained for Landsat 8 versus Landsat 7/8, Sentinel 2A/2B, and Terra/Aqua MODIS presented that the difference between these pairs was within 0.5–1% for most of the spectral bands. Total uncertainty obtained for these pairs of sensors using Monte Carlo Simulation varies from 2.5–4% for all bands except for SWIR2 bands, which vary up to 5%. The difference between EPICS Global and EPICS North Africa was calculated using the ratio of trend gain; the difference among them was within 0.5–1% difference on average for all the sensors and bands within a 0.5% uncertainty level difference

Calibração absoluta do radiômetro portátil CIMEL/CE313 em laboratório e avaliação das incertezas

Absolute calibration is required when the focus is the use of reliable quantitative data from sensors. This process determines the relation between the incident radiant flux and the digital number provided by the sensor. In this work we describe a methodology to perform laboratory absolute calibration of a sensor system, also presenting its uncertainty of measurement. The methodology was tested with the use of a portable radiometer CIMEL/CE313, without previous calibration information. This work represents a cross-calibration, involving two sensors: one of them used as reference. Thus, the radiometer CIMEL/CE313 was calibrated in relation to the ASD FieldSpec Pro spectroradiometer, recently calibrated by the manufacturer. The ASD FieldSpec Pro was considered the standard reference. Results showed that the relation of radiance values and digital numbers is linear. Based on this, we determined the radiometric calibration coefficients, gain and offset, which allow the conversion of digital number into radiance values. Moreover, the parameters offset, for bands 2, 3 and 4 of radiometer CIMEL/CE313, were consistent with zero value, when considering the uncertainties. So, only the gain parameter allows the conversion of digital number, generated by the radiometer into radiance values for these bands.Pages: 8962-896

Hyperspectral Empirical Absolute Calibration Model Using Libya 4 Pseudo Invariant Calibration Site

The objective of this paper is to find an empirical hyperspectral absolute calibration model using Libya 4 pseudo invariant calibration site (PICS). The approach involves using the Landsat 8 (L8) Operational Land Imager (OLI) as the reference radiometer and using Earth Observing One (EO-1) Hyperion, with a spectral resolution of 10 nm as a hyperspectral source. This model utilizes data from a region of interest (ROI) in an “optimal region” of 3% temporal, spatial, and spectral stability within the Libya 4 PICS. It uses an improved, simple, empirical, hyperspectral Bidirectional Reflectance Distribution function (BRDF) model accounting for four angles: solar zenith and azimuth, and view zenith and azimuth angles. This model can perform absolute calibration in 1 nm spectral resolution by predicting TOA reflectance in all existing spectral bands of the sensors. The resultant model was validated with image data acquired from satellite sensors such as Landsat 7, Sentinel 2A, and Sentinel 2B, Terra MODIS, Aqua MODIS, from their launch date to 2020. These satellite sensors differ in terms of the width of their spectral bandpass, overpass time, off-nadir viewing capabilities, spatial resolution, and temporal revisit time, etc. The result demonstrates the efficacy of the proposed model has an accuracy of the order of 3% with a precision of about 3% for the nadir viewing sensors (with view zenith angle up to 5°) used in the study. For the off-nadir viewing satellites with view zenith angle up to 20°, it can have an estimated accuracy of 6% and precision of 4%

The Development of Dark Hyperspectral Absolute Calibration Model Using Extended Pseudo Invariant Calibration Sites at a Global Scale: Dark EPICS-Global

This research aimed to develop a novel dark hyperspectral absolute calibration (DAHAC) model using stable dark targets of “Global Cluster-36” (GC-36), one of the clusters from the “300 Class Global Classification”. The stable dark sites were identified from GC-36 called “Dark EPICS-Global” covering the surface types viz. dark rock, volcanic area, and dark sand. The Dark EPICS-Global shows a temporal variation of 0.02 unit reflectance. This work used the Landsat-8 (L8) Operational Land Imager (OLI), Sentinel-2A (S2A) Multispectral Instrument (MSI), and Earth Observing One (EO-1) Hyperion data for the DAHAC model development, where well-calibrated L8 and S2A were used as the reference sensors, while EO-1 Hyperion with a 10 nm spectral resolution was used as a hyperspectral library. The dark hyperspectral dataset (DaHD) was generated by combining the normalized hyperspectral profile of L8 and S2A for the DAHAC model development. The DAHAC model developed in this study takes into account the solar zenith and azimuth angles, as well as the view zenith and azimuth angles in Cartesian coordinates form. This model is capable of predicting TOA reflectance in all existing spectral bands of any sensor. The DAHAC model was then validated with the Landsat-7 (L7), Landsat-9 (L9), and Sentinel-2B (S2B) satellites from their launch dates to March 2022. These satellite sensors vary in terms of their spectral resolution, equatorial crossing time, spatial resolution, etc. The comparison between the DAHAC model and satellite measurements showed an accuracy within 0.01 unit reflectance across the overall spectral band. The proposed DAHAC model uncertainty level was determined using Monte Carlo simulation and found to be 0.04 and 0.05 unit reflectance for the VNIR and SWIR channels, respectively. The DAHAC model double ratio was used as a tool to perform the inter-comparison between two satellites. The sensor inter-comparison results for L8 and L9 showed a 2% difference and 1% for S2A and S2B across all spectral bands

Method of Validating Satellite Surface Reflectance Product Using Empirical Line Method

Atmospherically corrected surface reflectance (SR) products are used for reliable monitoring of land surfaces and are the standard products of Landsat sensors. Due to increased demand for SR products, a need exists to verify that the L2C2 (Level-2 Collection-2) SR products are precise and accurate. The Level-2 Collection 2 (L2C2) SR Product has processed satellite imagery data that correct for atmospheric effects such as absorption and scattering, providing a more accurate representation of Earth’s surface. The validation of SR products using ground truth measurement is essential. This study aims to develop and evaluate a validation methodology for satellite SR products. Thus, the Empirical Line Method (ELM) is used here for atmospheric validation of remotely sensed data. Validation is performed using the SR derived from ELM tied to ground truth measurement. Absolute surface reflectance models of Algodones Dunes and the Salton Sea located in North America Sonoran Desert are developed to extend the temporally limited ground truth measurements. This model can give ground truth reflectance in any time frame independent of time constraints. The result of the absolute surface reflectance model of Algodones Dunes indicates that the model predicts the response of Algodones Dunes with an average accuracy of 0.0041 and precision of 0.0063 and gives ground measurements across all multispectral between 350 and 2500 nm. For the Salton Sea, the model predicts the response of the Salton Sea with mean absolute error (MAE) of 0.0035 and gives ground measurements across all multispectral between 350 and 2500 nm. The ELM generates atmospheric coefficients (gain and bias), which are applied to an image to obtain SR. Validation results indicated that for L9-OLI-2, L8-OLI, and L5-TM-SR products, the RMSE range is 0.0019 to 0.0106, 0.0019 to 0.0148 and 0.0026 to 0.0045 reflectance unit, respectively, and accuracy is within 0.0038, 0.0022, and 0.0055 reflectance unit across all spectral bands of L9, L8, and L5, respectively. On average, the validation result showed a strong linear relation between the L2C2 SR products and ELM SR within 0.5 to 1 reflectance units. These results demonstrate the high accuracy and reliability of the L2C2 SR product, providing valuable information for a wide range of remote sensing applications, including land cover and land use mapping, vegetation monitoring, and climate change studies